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Abstract 

We study the Moyal evolution of the canonical position and momentum variables. We compare it 
with the classical evolution and show that, contrary to what is commonly found in the literature, the 
two dynamics do not coincide. We prove that this divergence is quite general by studying Hamiltonians 
of the form p^ /2ra + V{q). Several alternative formulations of Moyal dynamics are then suggested. 
We introduce the concept of starfunction and use it to reformulate the Moyal equations in terms of 
a system of ordinary differential equations on the noncommutative Moyal plane. We then use this 
formulation to study the semiclassical expansion of Moyal trajectories, which is cast in terms of a (order 
by order in fi) recursive hierarchy of i) first order partial differential equations as well as ii) systems of 
first order ordinary differential equations. The latter formulation is derived independently for analytic 
Hamiltonians as well as for the more general case of smooth local integrable ones. We present various 
examples illustrating these results. 

1 Introduction 

In recent years there has been a lot of work devoted to noncommutative structures in mathematics and 
physics [1], [2]. In most cases one deals with the Moyal product, which first appeared in the context of 
deformation quantization [3, 4, 5, 6]. If we consider the position and momentum operators generators of 
the Heisenberg algebra A, then we can define the associated universal enveloping algebra in quantum phase 
space Ui, {A) , where c- functions are multiplied by resorting to a twisted convolution which admits a formal 
7i-expansion^ : 



f{Q,p)*g{Q,p) = exp 



ih / d d d d 
2 \dq dp' dp dq' 



f{q,p)g{q',p . 



(1) 

{,q\p')=iq,p) 
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^We shall only consider one-dimensional systems in this work. The generalization to higher dimensions is straightforward. 



or an alternative kernel representation valid in the set of phase space square integrable functions L2{C, dT), 
dT = dqdp [7]: 

f{q,p)*g{q,p) = j^jdq'dp'dq"dp" f{q',p')g{q",p")exp |-| [p{q' - q") +p'{q" - q) +p"{q - q')] 

(2) 
This is the so-called Moyal (or Groenewold) product. One can also construct a Lie algebraic bracket - the 
Moyal bracket - according to: 
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[f{q,p),9{q,p)] = ^sm 



h / d d d d 
2 \dq dp' dp dq' 



f{Q,p)9{Q',p' 



(3) 

(q',p') = {q,p) 



These provides the basic structure for the formulation of quantum mechanics in the phase space [7]- [18]. 
One of the remarkable properties of these expressions is the fact that they constitute formal deformations 
of the usual commutative product and of the Poisson bracket, respectively: 

f{Q,p)'^9{q,p) = f(.q,p) ■ 9(.q,p) + 0{h) 

(4) 
[f{q,p),9{q,p)] = {f{q,p),9{q,p)} + 0{n^) 

The quantum phase space endowed with the Moyal bracket admits a Heisenberg evolution for (real) c- 
functions according to: 

A{q,p,t) = [Aiq,p,t),Hiq,p)]^gj^^, Aiq,p,t = to) = Aoiq,p) (5) 

where H{q,p) is the Hamiltonian symbol (taken to be real and ?i-independent). In particular, for the 
fundamental variables {q,p), we have: 

Q{q,p,t) = [Q{q,p,t),H{q,p)]^g^p^, P{q,p,t) = [P{q,p,t),H{q,p)]^^^^^y (6) 

subject to the initial conditions: 

Q{q,p,t = to) =q, P{q,p,t = to) =p. (7) 

The subscript {q,p) in eq.(6) stresses the fact that the Moyal bracket is to be evaluated with respect to 
the canonical set {q,p). From (6,7) we may then write at time t = tQ-. 

Qiq,p,t = to) = -^Hiq,p), P(q,p,t = to) = -—H{q,p). (8) 

Equations (8) are reminiscent of the classical Hamilton equations: 

Qiq,P,t) = -^H{Q,P), P{q,p,t) = -^H{Q,P). (9) 

It is commonly claimed that, because of eqs.(8), the fundamental {q,p) evolve classically. The difference 
between Moyal and classical evolution would then only become apparent once functions of q's and p's are 
considered. 

The subject of Moyal dynamics is a subtle one and previous work in the field is scarce. In [9] the 
difference between classical and Moyal dynamics for general observables is clearly acknowledged, although 
no concrete statement is made concerning the canonical variables. In [19] we pointed out explicitly that, 
in the Moyal framework, Hamilton's equations for the fundamental variables (8) should be solved in the 



noncommutative algebra U^ [A) (in which case they become vahd at later times) . The formulation of 
Weyl-Wigner quantum mechanics in the Heisenberg picture (a key ingredient for the discussion of Moyal 
trajectories) has been studied in many papers (see for instance [14, 20]). One of the most significative 
papers in the field of Moyal dynamics is [20], where the Moyal correction to the classical Hamiltonian 
equations is derived order by order in fi by using a cluster-graph representation of the star-exponential 
as well as a recursive hierarchy of classical transport equations (see also [21, 22]). We shall review this 
method in section 3.1. 

The purpose of this paper is to further clarify certain aspects of Moyal dynamics and present some new 
results. In particular we shall prove with concrete examples that the claim about the coincidence of Moyal 
and classical trajectories is false for most Hamiltonians. The subject of trajectories generated by a bracket 
on a symplectic (or more generally a Poisson manifold) with a star product is a notoriously difficult one. The 
present work is devoted to the simpler case of the flat phase space of a one dimensional system T*M ~ IR 
with a global Darboux chart and symplectic form oj = dqf\dp. We shall elaborate on the specific features of 
Moyal trajectories which make them distinct from their classical (Hamilton) counterparts, namely the facts 
that, in general, Moyal evolution does not act as a coordinate transformation nor as a symplectomorphism 
in phase space. Moreover we propose strategies for solving the equations, if not completely then at least 
formally or semiclassically (in a fi power series). These results are developed using noncommutative phase 
space methods while never resorting to operator methods. 

We feel that Moyal dynamics is an interesting and challenging mathematical problem per se. However 
Moyal trajectories are also expected to play an important role in certain physical situations. In ref.[23] 
we proved that the expectation value of an observable A (typically the position operator) in a state with 
density matrix p is given by: 



A{t) = Tr {A{t)p) = F 



-^-q;] ^M{z,t) 



(10) 



In the previous equation z = {q,p), -^ = (^)^)) AM{z,t) is the Moyal evolution of A{z) (in other 

words, the Weyl symbol of the Heisenberg operator A{t) ) which is the solution of eq.(5) and F{a) is the 
symplectic Fourier transform (or chord function) of the Wigner function F{z): 



F{a) 



dz F(z)e 



icr{a,z) 



a[a, z) 



MP 



If the state of the system is described by a coherent state with wavefunction 



:ii) 



i^aix) =[ — ] exp 



mu) 



T^i^-Q) 



2h 



H [ X 

n \ 2 



a 



mijj ip 



then eq.(lO) reduces to: 



A{t) = exp 



fi d^ hmuj d'^ 



Amoj dq^ 



+ 



4 dp^ 



^M{q,P,t). 



(12) 



(13) 



This expression is particularly well suited for the semiclassical expansions in powers of Planck's constant. 
It is well known that expectation values such as (13) (or more generally (10)) play an important role in 
quantum optics [8] and in determining effective actions in quantum mechanics and quantum field theories 
[24], [25]. Moreover the method of trajectories finds important applications in quantum hydrodynamics 
[26] and transport models in quantum chemistry and heavy-ion collisions [27] and it is also expected to be 
crucial if one aims at fully grasping the nonlocal nature of noncommutative quantum field theories [1]. 

We shall start by recapitulating some well known features of classical dynamics, namely that classical 
evolution is a canonical and a coordinate transformation, and that the Hamiltonian is a constant of motion 



when not explicitly time dependent [28], [29] (section 2). We shall then show how these facts translate to 
Moyal dynamics (section 2) and propose new formulations of Moyal equations and new strategies to solve 
them (section 3). All the results will be illustrated with examples (section 4). The Hamiltonian of the 
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form H = -^ + V{q) will be studied in detail (section 5). 

2 Classical and Moyal dynamics 

Some of the results of this section are well known. We present them here for completeness and to highlight 
the differences between Moyal and classical evolutions. We will present simple proofs of the results whenever 
we feel that there is some subtlety that should be emphasized. We shall use the subscripts C and M for 
classical and Moyal evolutions, respectively. 

One possible formulation of classical dynamics is in terms of a classical transport equation. Let then 
Ac{q,p,t) be the solution of the first order partial differential equation 

Ac{q,P,t) = {Ac{q,p,t),H{q,p)}(^g^p^ , (14) 

subject to the initial condition 

Ac{q,p,t = to)=Ao{q,p). (15) 

Then it is easy to prove that: 

Theorem 2.1: The solution is 

Ac{q,p,t) = Ao{Qc{q,p,t),Pc{q,p,t)) , (16) 

where {Qc{q,p,t), Pc{q,p,t)) are the solutions of Hamilton's equations (9) with initial conditions (7). 
By resorting to Jacobi's identity and theorem 2.1, we can also show that: 

Theorem 2.2: The transformation 

{q,p) ^{Q = Qc{q,P,t),P = Pc{q,P,t)) (17) 

is canonical, i.e. 

{A{Qc{q,p,t),Pc{q,p,t)) ,B {Qc{q,p,t),Pc{q,p,t))}(^g^p^ 
= {A{Qc{q,p, t), Pc{q,P, t)),B{Qc{q,p, t),Pc{q,p, t))}(Qe,Pc) ' (^8) 

for general phase space functions A, B. 
Lemma 2.3: The Hamiltonian is a constant of motion. 
Proof: For an arbitrary observable we have: 

Ac{q,P,t) = {Ac{q,p,t),Hoiq,p)}^g^p) = F{q,p,t), (19) 

where we introduced the notation HQ{q,p) = H{q,p,t = tg)- Differentiating the previous equation with 
respect to time, we obtain: 

F{q,P,t) = {AcHo} + {AcHo} = {Fiq,p,t),Ho{q,p)}^g^j,y (20) 

From theorem 2.1 it then follows: 

F{q,p, t)=F {Qc{q,P, t), Pc{q,P, t),t = to) . (21) 



In particular, let us choose A = H. Then, since F{q,p,t = to) = 0, we conclude from the previous analysis 
that F{q,p,t) = and from (19): 

Hiq,p,t) = {Hiq,p,t),Hoiq,p)}(^g^p) = O.a (22) 

An obvious corollary of the previous result is the following: 

Corollary 2.4: The functional form of the Hamiltonian remains unchanged in the course of classical 
evolution: 

H {Qc{q,p,t),Pc{q,p,t)) = H{q,p). (23) 

The next theorem is the crux of the matter when we compare classical and Moyal evolutions. It states 
that the equations of motion are the same at any time t: 

Theorem 2.5: The partial differential equations (14) for the position and momentum are equivalent to 
the set of ordinary differential equations: 

a IT f)ff 

Q = {Q, H}^Q,p) = ^, P = {P, H}^Q^^^ = - — . (24) 

Proof: We prove the statement for Q. From eq.(14), theorem 2.2, lemma 2.3 and its corollary we get: 

Q{q,P,t) = {Q{q,p,t),H{q,p)}^g^^ = {Q{q,p,t),H {Q{q,p,t),P{q,p,t))}^g^^ = 

(25) 
= {Q{q,P,t),H{Q{q,p,t),P{q,p,t))}^QP^ ^^ Q = {Q , H {Q , P)} ^^q p^ , 

where we omitted the subscript " C" for simplicity and used eq.(18).n 

Remark 2.6: Notice that this presentation of classical dynamics is somewhat unusual. We posited the 
evolution to be given by the classical transport equation equation (14) even if Ac stands for the position 
or the momentum. For instance: 

Qc{q,P,t) = {Qc{q,P,t),H{q,p)} = -^— q~~q~' Qc{q,P,t = to) = q, (26) 

and a similar equation for Pc{q,p,t). This is a partial differential equation. It is only through theorem 
2.5 that we show the coincidence of the solutions of the partial differential eq.(26) and of the system (24) 
of ordinary differential equations. We have deliberately chosen this presentation to render the passage to 
Moyal dynamics more transparent. As a word of caution we stress the double role played by the position 
and momentum in this approach. They are simultaneously the variables with respect to which one performs 
the partial derivatives in say eq.(14) as well as the initial conditions (7). 

Let us now turn to Moyal dynamics (6,7). It is well known that 

Theorem 2.7: The formal solution of the Moyal equation is given by: 

AM{q,p,t) = U{t)*Ao{q,p)*U{-t), U{t) ^ ef-*°)^('''P)/^ (27) 

Here 4>{l3,q,p) = e* ^'^ is the noncommutative exponential, solution of: 

^ = B^(t) = (l)^B, </,(/3 = 0,g,p) = l, V((z,p)Gr*M~IR2. (28) 



Remark 2.8: Contrary to what happens in classical mechanics (cf.eq.(16)), Moyal evolution does not act 
as a coordinate (local) transformation. In general: 

AM{q,P,t) ^ Ao{QM{q,P,t),PM{q,P,t)) . (29) 

The exception are flows generated by Hamiltonians which are at most polynomials of degree 2 in {q,p). 
We shall construct an explicit example which corroborates (29) in due course (see remark 4.4 below). 
An immediate consequence of theorem 2.7 is the analog of Lemma 2.3: 

Lemma 2.9: The Hamiltonian is conserved under Moyal evolution. 

The following, well known result, is the Moyal counterpart of theorem 2.2: 

Theorem 2.10: Moyal evolution 

iq,p) ^ (Q = QMiq,P,t),P = PMiq,P,t)) (30) 

being a unitary transformation, preserves the Moyal bracket: 

[QM{q,P,t),PM{q,P,t)](^q^p) = 1. (31) 

Conversely, all solutions of the constraint (31) are unitary transformations [7]. 

Theorem 2.10 states that it is the Moyal bracket rather than the Poisson bracket that is preserved under 
time evolution. Let us clarify this point. There are, obviously, transformations which are both canonical 
and unitary. 

Example 2.11: A transformation which is both canonical and unitary. Let us consider the symplectic 
transformation: 

Q{q,p) =aq + bp, P{q,p) = cq + dp, (32) 

where: 

de.(°;;)=l. (33) 

This transformation is both canonical as well as unitary: 

[Q'^]fep) = W'^}fep) = l- (34) 

However, there are transformations which are unitary but not canonical and vice-versa. 

Example 2.12: A Transformation which is unitary but not canonical. Let us deflne the transformation: 



Q(g,p)=/3e«/^ P{q,p) = e-'^'^ 



p + 7 sinh f — - — j 



(35) 



where /3, 7 are positive real constants with dimensions of length and momentum, respectively. This 
transformation is invertible by the inverse function theorem. The Jacobian reads: 

d{Q,P) .^ p. ,^ 2/377r (2(57:p\ 

-Wp^ = ^^' ^^(-) = ' + ^ ^°^^ [^) > ' ('') 

We conclude that this is not a canonical transformation. On the other hand we have: 

d d d d\ r . n /2^\2"+i /2/37rp\ 

TTTT ~ TTTT] P{Q^P) = PI {^-] ^osh ^ ^^ , n>l 



'^^^•"^^o-.Tp-o-pir, p(^-p) = i^-'[ir) ™Hirj' "2'- p^' 



where the derivatives d and d act on Q and P, respectively. Consequently: 

2 / ^ ~* ^ ^ \ 2n+l 

= 1 H — Y' cosh ( ^^ j sin tt = 1. 



(38) 



And so this is a unitary transformation. 

As the latter example shows, unitary transformations are not canonical in general. Likewise, canonical 
transformations are not unitary in general (see example 1 below). Therefore, one should expect Moyal and 
classical trajectories to diverge. Indeed, this is stated in the following theorem. 

Theorem 2.13: At an arbitrary instant t ^ to we have, in general, for Q = Qm , P = Pm'- 

Q{t) / [Q(t),i7(Q,P)](Qp) = -QpH{Q,P), P{t) / [P(t),i7(g,P)](Qp) = -—HiQ,P) (39) 

Proof: For Ao{q,p) = q the solution of eq.(5) takes the form (cf.(27)) Q{q,p,t) = U{t) -kq-k U{—t). And 
so: 

Q{q,p, t) = U{t) * [q, H{q,p)\^^^^ * U{-t) = U{t) • F,{q,p) * U{-t), (40) 

where Fq{q,p) = dH/dp. Unless Fq{q,p) is linear, we have (cf. remark 2.8): 

U{t) * Fg{q,p) * U{-t) / Fg {Q{q,p, t), P{q,p, t)) . (41) 

Thus, in general: 

Q(.t)^^H{Q,P). (42) 

A similar argument holds for PM{q,p,t)-n 

3 Other approaches to Moyal dynamics 

The main conclusion of the previous section is that the partial differential equations (6) governing Moyal 
dynamics are not equivalent to the system of ordinary differential equations (9). Since eqs.(6) are of 
arbitrary order (possibly infinite) the problem of solving them is notoriously difficult. In this section 
we investigate various alternative formulations of Moyal mechanics. One obvious approach is that of 
evaluating semiclassical corrections (order by order in h) to the classical solutions. In section 3.1 Moyal 
dynamics is written in terms of a ?i-hierarchy of recursive first order linear partial differential equations. 
This hierarchy has been previously presented in [20]. Here we re-derive it in other (possibly simpler) terms. 
In section 3.2 Moyal equations are re-written as a system of ordinary differential equations in the space 
of noncommutative phase-space functions. This is an interesting formulation allowing, in some cases, to 
solve Moyal equations exactly (see section 4). Furthermore, and as a by product, Moyal dynamics will be 
written in terms of a /i-hierarchy of systems of ordinary first order linear and inhomogeneous differential 
equations in section 3.3 and 3.4. In a first step this result is derived for analitic Hamiltonians (section 3.3) 
and then generalized to the case of smooth, locally integrable ones (section 3.4). In this section we shall 
write Qm = Q, Pm = P for simplicity. 



3.1 A hierarchy of first order partial differential equations 

If we regard Moyal evolution as a partial differential equation, then we can easily derive an equivalent 
infinite hierarchy of linear inhomogeneous partial differential equations. Indeed let Q{q,p,t) and P{q,p,t) 
be the solutions of (6,7). A key feature of these solutions is that they display a regular behavior on h and 
admit a formal asymptotic series: 



+00 

E 

n=0 



Q{q,P,t) = J2 ?i"Q^"H9,P,*), P{(l,P,t) = E n''P^''HQ,P,t), 



(43) 



ra=0 



where, obviously, Q^^'{q,p,t), P^^'{q,p,t) are the classical solutions. A simple analysis reveals that only 
the even order terms are non-vanishing, provided the Hamiltonian is independent of h. If we substitute 
(43) in the Moyal equations, we get the classical Hamilton equations for n = and: 



p(2n) =^^^Jp(2fc),i7' 






2{n-k) 
2(n-fc) 



n> 1 



n > 1 



(44) 



where we used the notation: 
[A{q,p),B{q,p)]^^ 



-I) 



d d d d 



2ra+l 



M<l,P)B{q' ,p' 



(2n + l)!4" \dqdp' dp dq' 

Here Q^'^^\t) and P^'^'"-\t) satisfy the initial conditions: 

Q(2")(g,p,t = to) = P^^''Hq,P,t = to) = 0, Vn > 1, V(g,p) G IR2. 



{q',p') = iq,p) 



(45) 



(46) 



The important thing to remark in eqs.(44) is the fact that they constitute an infinite hierarchy of linear 
inhomogeneous partial differential equations. For each order 2n they are in fact a system of two independent 
Hamiltonian transport equations with inhomogeneous terms: 



n— 1 



n—1 



(47) 



(48) 



Ft'HQ,P,t) = E[Q^''^^]2(™-fc) and F('-Hq,p,t) = Y.[P^''\ H],^^^,) 
fc=0 fc=0 

exclusively dependent on the solutions of order 0, 2, • • • , 2n — 2. The solution of eq.(47) is then [20]: 

Q(^^)iq,p,t) = JldTFt\Q^'Hq,P,r),P^^Hq,p,T),T), 

p('^^\q,p,t)=JldTFi'''\Q(^\q,p,T),P^^Hq,p,T),T). 

3.2 Ordinary differential equations in the non-commutative phase space 

We start by introducing the concept of •-function. Intuitively, these are functions whose expressions 
are written in terms of ^-products, i.e. the fundamental operations used to express the •-functions are 
summations, multiplication by scalars and •-products. The most interesting feature of these functions is 
that their functional form is preserved through Moyal time evolution. Hence they provide the right objects 
to cast eqs.(8) in an invariant form. 



The main question is then what kind of functions / : T*M — > C can be written as ^-functions. To 
begin with the simplest case, let us consider an analytic function N{q,p) in the phase space variables. It 
is trivial to check that [19] (let Oi = q and O2 = p)'- 



°° 1 9"A^ 

n=0 l<ii,...,j„<2 '1 ' 



iO,,*...*0,Js^Wi<l^P)]s (49) 

(g,p)=(0,0) 



where S stands for total symmetrization. We used the fact that (q-k ... -kq-kp-k ... •kp)s = q^p^ and 

^ ^ ^ ^ ^ ^ 

n terms m terms 
introduced the notation [N{q,p)]s to denote the •-function associated to N{q,p) which is written in a 

completely symmetric form. Notice that the objects [N{q,p)]s are naturally characterized as elements of 
the enveloping algebra U-^ (A) of the Heisenberg-Weyl Lie algebra, i.e. as series in the non-commutative 
coordinates. 

Let us now extend the concept of •-function beyond the case of phase-space analytic functions. We 
start by noticing that an arbitrary locally integrable function F{q,p) admits a (generalized) Fourier repre- 
sentation: 

F{q,p) = fd^fdri a(e,r/)e^««+^'?P. (50) 



in terms of possibly distributional coefficients a(^,r/). The set of functions of the form (50) comprises the 
set of all phase space analytic functions as well as the set of square integrable functions (in fact formula 
(50) provides a representation for generalized functions). 

All functions of the form (50) can now be written as •-functions: 



[F{q,p)]s = jdCjdV a{C,7^)el^i+^^P. (51) 

For example, if F{q,p) = q'^p, then we have a{S,,T]) = —i6"{S,)6'{7]). It then follows that [F{q,p)]g = 
{q-kq-kp + qkcp-kq + pkcqkcq)/3. Notice that the identity F{q, p) = [F{q, p)]s follows immediately from: 

^iS,q+irip _ ^i^q ^ ^ivP^^S.V — g«C9+«^P^ C52) 

where we used the Baker-Campbell -Hausdorff (BCH) formula. 

To proceed we notice that •-functions can be cast in many different functional forms which is a direct 
consequence of ordering ambiguities (e.g. q-k p = p-k q + ifi). An important alternative to the complete 
symmetrization is the standard-antistandard symmetrization that we shall use in section 4 to study a 
specific example. In this case we write a combination of a term, where all the g's stand to the left and its 
conjugate. For F{q,p) = q^p, we have [F{(l-,py\sAS ~ {Q'^ Q'^P + P'^ Q'^ o) l"^- This can be expressed as: 

[F{q,p)]sAS = lJd^Jdv a(e, v) sec (|^) {e^' * e^ + e^ * e^^^) . (53) 

Lemma 3.1: For a phase space function satisfying (50), the following holds: 

[F{q,p)]sAS = F{q,p). (54) 

Proof: Again using the BCH formula, we have: 

^JdCjdr] a(e, v) {el^" * eT + el'"' • e^^") = ^ J d^' di] a(^, r/)el«''+*''P (^e^^Hv/^ + eiH^/2^^ ^ 

(55) 
= Jdadr] a(e,r/)e^«^+*'?Pcos(^) =cos(|g|^)F((7,p) 



from where the identity (53) follows. n 

Now let us go back to the equations of motion. At time t = tg eqs.(8) hold. However, they are only 
valid at t = to- Let U{t) be the symbol of the evolution operator and let [Fq{q,p)]^ be the ^-function 
associated to Fq{q,p) = ^ and written in the order O = S* or 5^5". We then have: 

Q{t) = [Q{t),H](^g^p)=U{t)^[Q{0),H](^g^p)*U{-t) 

= Uit)*[F,iq,p)]^*U{-t) = [F,{Q{t),Pmo^ (56) 

and a similar equation for the momentum. The important thing to remark is the fact that the equations of 
motion (56) are now valid for all times in some interval ]tQ — €,tQ + €[ (e > 0). We may now regard these 
equations as a system of ordinary differential equations for the position and momentum. We just have to 
bear in mind that we have to look for solutions in the set of *- functions. The example of section 4 will 
illustrate this procedure and show how, in some cases, the formulation (56) allows for the exact solution 
of Moyal dynamics. . 

3.3 A hierarchy of ordinary differential equations for analytic Hamiltonians 

As a direct application of equation (56) we now study the 7l-semiclassical expansion of the Moyal evolution 
of position and momentum (43). Let us consider the case where the Hamiltonian H{q,p) is a real analytic 
function, independent of h. In this case [Fq{q,p)]s can be written as: 

oo 

[Fgiq,p)]s= E C^n,™(g>"*') (57) 

n,m=0 

where k = n + m — 1 for (n, m) ^ (0, 0) and k = for n = m = 0, Cn,m are the Taylor coefficients and: 

{qy"^^k)^{q^...^q^p^...^p)s. (58) 

n terms m terms 
Using this notation eq.(56) is written: 

oo 

g= E Cn,m(Q"^™*') (59) 

n,m=0 

Let us then expand this equation in powers of h. We consider the expansions (43) and introduce the 
additional notation: 



^ 1 /7\"//9 8 8 8 \"- 



n=0 

It then follows that: 
where 



• (60) 

(q' ,p') = {q,p) 



(^qnpm^k^ = E fl'\Q''P"'i^% (61) 

r=0 



(62) 



2_/a=l "ia + 2_//3=l Jl3 + Z_/7=l '7 — T 
ia,Jl3,ly G {0, ...,r} 
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s 



and so the leading order eq.(59,62) is the classical Hamiltonian equation, and to order 0{h^), r > 1 it 
reads: 

oo 

qm = j2 a,m(g"p"*')r (63) 

n,m=0 

oo 

n,m=0 

together with the initial conditions: 

QM(t = to) = P^^)(t = to) = 0, r>l (64) 

and the inhomogeneous term is: 

/,(QW,...,Q(^-1),PW,...,P(^-1)) (65) 

oo 
= E Cn,m E {Q^''^ *h Q^''^ */.••• *l.-. Q^'"^ ^Ir. P^''^ *^„+i •••*/. P^'-^ 

ia,jf3 G {0, ...,r - 1};/^ G {0, ...r} 

A similar set of equations are valid for P. Finally, it is easy to verify that the initial conditions will impose 
the odd contributions to the semiclassical expansions of Q, P to be identically zero. 

3.4 The hierarchy of ordinary differential equations in the general case 

We now generalize the previous approach to local integrable Hamiltonians (still real and time, h indepen- 
dent). 

Theorem 3.2: Let the Hamiltonian admit a Fourier representation as in (50). Then in the space of 
formal series [ [ ?i ] ] the Moyal evolution of position and momentum is dictated by a hierarchy of sets of 
inhomogeneous first order linear differential equations. 

Proof: Let us write z = {q,p), Z = (Q, P). The Moyal equations at time t = to read: 

Z,{t = to) = Fi, F, = J^jd,H{z), (66) 

where sum over repeated indices is understood, H is the Hamiltonian and we used the compact notation 

^* ~ (^' ^) ' "^11 ~ "^PP ~ ^' '^<1P ~ ~'^Pi ~ -'■• Under the assumptions of the theorem the vector function 
F admits the Fourier transform: 

F{z)= (di /(e)e^«•^ (67) 



where ^ = (^g, ^p) and ^- z = q^q + pS,p- Following the procedure discussed in section 3.2, we may write at 
time t > to: 

Z{t)= fd^ fiOef^^'^. (68) 



In ref.[15] we proved that the noncommutative exponential may be expanded in the form: 

ef = (l + h^A2 + h^A^ + • • •) e^, (69) 

where A2, A4, ■ ■ ■ depend only on B{z) and its derivatives. For instance: 



^2(2) = -JikJji (didjB) 



^(dkdiB) + ^idkB){diB) 
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(70) 



It should be noted that B may also depend upon h. If we expand the Moyal trajectories 

z{t) = z^'^^t) + n2z(2)(t) + n^z^^\t) + ■■■ 

as well as the coefficients of the exponential 



(71) 



A^,{Z) = 4J(^^°^) + ?i'4?(^^°U(')) + ^^4fc (Z(0), Z(2), Z(4)) + • • • , fc = 1, 2, ... (72) 

we get from (68): 

^(0)(t) + ;j2^(2)(t) + h^Z^^){t) + h^Z^^\t) + ■■■ = ! di /(O {e^«-(^(°'+?i^^(^>+?i^^(^'+?i^^(«>+-) 



+ 



+/i2 



.2.(2) 



4'^^(ZW) + ?i24^)(z(o), Z(2)) + n44"^(z(o), Z(2), zW) + 



4 4(4), 



„iC-(Z(0)+a2Z(2)+a2z(4) + ...) 



+ 



+n^ 



Af (ZW) + /i242)(z(0), Z(2)) + . . .1 Siz^''>+r^zi^)+-) ^ ^6 M0)(z(0)) 



+ ••• 



(73) 



And so to order 2n: 



aR 



^^Ho = / di friOe'^-'"' Oe • ^(^") + •••) = 4'"^(0 • ^(^^°H0) + 



dZ} 



(74) 



where the dots stand for terms which depend only on Z^^\ Z"-*, • • • , Z'^" '^'. As advertised these constitute 
an infinite hierarchy of ordinary first order linear differential equations. n 

Remark 3.3: The theorem is stated in the space of formal series, which means that questions of conver- 
gence are ignored (this was already the case in sections 3.1 and 3.3). The important thing to stress is 
the fact that the hierarchy of partial differential equations of section 3.1 has now been transformed into a 
hierarchy of ordinary differential equations. The proof also indicates how one can write the inhomogeneous 
term in the system of ordinary differential equations to each order 0{h "). For instance, to order 0{h ) 
the inhomogeneous term reads: 

Jd^ /(e)e^«-^'"Uf(Z(0)) = 

= -Idc f{Oe'^■^'"'J^kJJl[^^^,{ic■z(^))] X [3^afea^(i^z(o)) + 2^afc(i^z(o))9^(^^z(o)) 



-JikJjl (didjZa I X 



JQ [OkOlZ^i^ 



dk-diZ, 



^'^^ -^ + h (^.^f ^ (d^zl'^ 



dZ, 



>l^c 



92 



5^f^ 






Altogether we may then write: 



zP{t) = zf )(t)^(z(o)) - ±j.,j,, (5.a,zr)) {dM'^) ^^{z^'^) 
-y.,j,, (5.a,z(°)) (d,zf^) (d,zf^) ^^M^r^izio)). 

This equation could be equally derived using the method of section 3.3 for analytic Hamiltonians. 



(75) 



(76) 
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4 Example 1 

The following example shows explicitly that Moyal and classical evolutions are different. Moreover it will 
illustrate the formalism developed in section 3.2. We shall assume to = for simplicity and we start by 
computing the classical evolution generated by the following Hamiltonian with position dependent mass: 



H{q,p) 



q2p2 
imp' 



where / is some real constant with dimensions of length. The Hamilton equations read: 



Qc 



QcPc i, _ QcPh 



2mP 



Pc 



2ml 



2 ■ 



We conclude that Qc ■ Pc is a constant of motion: 

Qc{q,p,t) ■ Pc{q,P,t) =q-p, 
Substituting this expression in the first equation (78) we get: 



Vt. 



Qc{t) 



qp 

2mP 



Qc{t). 



And so, from (7): 



Upon substitution in (79), we get: 



Qc{q,P,t) = qe'^ 



Pciq,P,t) =pe 



IP* 



ipt 



VtG K. 



VtG R. 



(77) 

(78) 

(79) 

(80) 
(81) 

(82) 



By construction this is a canonical transformation at all times. However, this is not a unitary transforma- 
tion. To check this, let us first compute Qc * Pc using the kernel representation (2). If we Wick rotate 
t -^ ir, we get after some integrations and an integration by parts: 



Qc(.q,P,t)*(q,p) Pc{q,P,t) 



qp + ih/2 



q *(q,p) P 



, , f ht V 



n2 



l+(^i_^ 



TinP J 



(83) 



It then follows that: 

[Qc{q,p,t),Pc{q,p,t)] 



M - n' 



I'm[Q*(q,p)P 



1 + 



/ ht 



\4mf 



/I, Vt / 0. 



(84) 



And so this is not a unitary transformation. From (31) we conclude that the Moyal evolution must differ 
from the classical solution. The Moyal evolution for an arbitrary observable A is dictated by the equation: 



A 



M 



qp 



mL 



d d\ h 

qiT -P^r + 



^3 



P-. 



Q3 



P \ dq dp J 8mP \ dq^dp dqdp 



A 



M- 



(85) 



Now this can be solved for {Qm,Pm) with the initial conditions (7). This is feasible in this case. But, in 
general, it is the fact that we are dealing with a partial differential equation that makes Moyal evolution 
difficult to solve. We shall follow the alternative route presented in section 3.2. This approach is reminiscent 
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of the classical solution presented above. We shall then solve an ordinary differential equation in the 
noncommutative space of •-functions. The equations of motion at time t = read (of. (8)): 

These equations are valid only at time t = 0. To write the corresponding ordinary differential equations 
at another time t we will have to write the right-hand sides in terms of •-products. We resort to lemma 
3.1 and write the equations (86) in the SAS order: 

Q = j:^{Q*q*p + p*q*q), p = -j:^{Q*p*p + p*p*q)- (87) 

If U{t) is as in (27), then if we act with U{t)-k on the left and with -kU{—t) on the right of these equations, 
we get: 



'3m = -; — ^ {Qm * Qm * Pm + Pm * Qm * Qm) , Pm = —-. — Y) {Qm * Pm * Pm + Pm * Pm * Q 



M 



Amr 4m/ 

The important thing to remark is that, contrary to eqs.(87), the previous expressions hold for all t G ] — e, e [ 
with e > 0^. We may thus regard eqs.(88) as ordinary differential equations for {Qm, Pm)- We just have to 
bear in mind the fact that we have to find solutions in the space of •-functions . Since [QM{t), Pnit)] = 1 
at any time t, we have: 

3j {Qm * Pm) = Qm * Pm + Qm * Pm = -^^^ {Pm * Qm * Qm * Pm — Qm * Pm * Pm * Qm) = 

(89) 
= 5;^ {Pm * Qm * Qm * Pm — Pm * Qm * Qm * Pm + ifiPM * Qm — ihPM * Qm) = 0. 

And so, as in the classical case, Qm * Pm is a constant of motion: 

III 
QM{q,P,t) •(g^p) PM{q,P,t) =q-kp = qp+—. (90) 



Substituting (90) in the first eq. 



1 

4mP 



QM{t) = 7372 \.QM{t) * {qp) + {qp) • QM{t)] ■ (91) 



The solution of this equation subject to the initial condition (7) is: 

tqp tqp 

QM{q,P,t) = eJ^ *q^er^. (92) 

Similarly, we get for Pm'- 

PM{q,P,t) = e^^^ *p*e^^^. (93) 

By construction this is a unitary transformation. To finish our calculation, we just need to compute 
explicitly the •-products in (92,93). The noncomutative exponentials have been computed for quadratic 
functions in [15], [9]. The result is: 



, , , ±-^ ^f fit \ 

cP±{t,q,p)^e."-' =sec^^^^j 



exp 



2 / ht 



ziz—qptan 



h KSml 



n 



2 



t| < ^^. (94) 



We show below that e ; 
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And so: 



QMiq,P,t) = (l)+-kq*(j)+= l+itan(g^j (^+ • (#+) 



sec'U^Uexp 



i'?Ptan(^) 



t < 



2-Kmf 

h ' 



(95) 



where the last ^-product was again evaluated using the kernel representation (2). A similar calculation 
leads to: 



PM{q,P,t) = sec (-— -^jpexp 



ht 



4mP 



2 / ht 



t < 



l-KmV 



h ' 



(96) 



Remark 4.1: The mapping {q,p) -^ {QM{t)^PM{t)) for this system constitutes another example of a 
transformation which is unitary but not canonical. Indeed we have: 

{QM(g,P,t),PM(g,P,t)}(,,rt = sec' [-^;^) /I' V* e ] - ^— , ^— [ \{0} . (97) 



Remark 4.2: As expected in the limit h ^ 0, the previous solutions coincide with the classical ones 
{Qc{Q:P:t)iPc{q,P:t)) (cf.(81,82)). Morcover, we have checked explicitly that the solutions to order h 
Q^'^\t),P^'^'{t)] satisfy the ordinary differential equations (76). 



Remark 4.3: Finding solutions in the noncommutative space of •-functions might be more restrictive 
than in classical mechanics. In our example, the Moyal solutions are only valid on the bounded interval 
t G ] — ^2™_ ^ 27rmi_ J ^ whereas the classical solutions are valid everywhere on [R . This is a very important 
point. It highlights the danger of relying too heavily upon /i-expansions. Indeed, if we consider say the 
second order solution for the position. 



Quii^P^'t) = Qc{q,p,t)- 



1 + ^ 

6mP J 



we realize that it is a smooth function everywhere on T*M ~ IR and for all t G IR. The same can be 
said about all the remaining orders. However, our exact solution (95) is singular at |t| = ^™ in sharp 
contrast with the perturbative result. 

Remark 4.4: This example illustrates the claim that Moyal evolution does not act as a coordinate (local) 
transformation. Let us consider the phase space variable: 



Ao{q,p) =q-kp= qp + 



ih 



(98) 



Since this is a constant of motion, we have: 



AM{q,P,t) = Uit)*Ao*U{-t) = QM{t)*PM{t) = qp + 



ih 



(99) 



Inverting (95,96) we get: 



q{QM,PM,t) = cos^ 



ht 

4mP 



Quexp -\QmPm sm 



2^ 1 ^os' 



ht 

4mT^ 



P (Qm, PM,t) = cos2 (^) Pm exp [^QmPm sin [^) cos^ (^) 



(100) 
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Upon substitution in (99) we obtain: 



Am = cos^ (^^) QmPm + y . (101) 

If Moyal evolution were a coordinate transformation we ought to have (cf(16)): 

AM{t) = Aq {QM(.t),PM{t)) = QmPm + y. (102) 

Notice also that eq.(lOO) is not the time reversal of eqs. (95,96) (i.e. cannot be obtained by making t -^ —t) 
as it would be the case if Moyal evolution were a coordinate transformation. 

5 Example 2 

The previous example is admittedly somewhat contrived. It has the virtue of being exactly solvable and of 
yielding distinct classical and Moyal evolutions. However, the claim that these evolutions are different in 
general remains valid even if we consider more realistic systems. For Hamiltonians which are polynomial 
of g's and p's of degree higher than 2, there will generally appear discrepancies between the classical and 
Moyal evolutions of order 0{h ). Indeed, let us consider a Hamiltonian of the form: 

2 

H{q,p) = ^ + Viq), (103) 

and define {z = {q,p)): 

' nf{z) ^ [z, H{z)]^ , ni{q,p) EE [ni{z),H{z)], , • • • 

_ kl{z)^{z,H{z)}^, Al{z)^{Kl{z),H{z)}^r-- 
A simple calculation reveals that: 



(104) 



nj{z)=Kl{zl . = l,---,5, f7^(^) = A^(^)-^F(3)(5)y(4)(g) (105) 

where V^^'{q) is the n-th derivative of V{q). Likewise for the momentum we get: 

n^^{z)=K^{z), ^ = l,---,4, al{z)=Kl{z)-^^V^^\q)V^'\q). (106) 

The solutions of the classical {Zc{z,t)) and of the Moyal {ZM{z,t)) equations of motion admit the formal 
expansions: 

oo , oo , 

Zciz, t) = z+J2 -,K{^)^ Zm{z, t)=z+J2 -,Ki^)- (107) 



_, n! ' — ' n\ 

n=l n=l 

By virtue of (105-107), we conclude that Qm diverges from Qc at order 0{t^): 

QMiq,P,t) = Qciq,P,t) - J^V^^\q)V^^\q) + 0{t'), (108) 

and Pm diverges from Pq at order 0{t^): 

PM{q,P,t) = Pc{q,P,t) - -^^V^^\q)V^'^\q) + ©(t^). (109) 
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For the quartic anharmonic oscillator V{q) = ■^muP'q^ + -^q^ (A > 0) and thus the previous equations yield: 

QM{q,P,t)=Qc{q,P,t)-——^q + 0{f), PM{q,P,t) = Pc{q,p,t) - ^ q + 0{t''). (110) 

It may appear from this analysis that a cubic potential generates a classical trajectory in Moyal mechanics. 
However, even this is not true. If we study the momentum to order 0{t'^) we obtain a term of the form: 

This example illustrates explicitly that the Moyal and classical evolutions of the fundamental variables 
{q,p) are generally different. 

6 Conclusions and outlook 

We have shown explicitly with various counter examples that, in general, Moyal and classical dynamics 
differ for the canonical variables of position and momentum. We argued that this stems from the fact 
that Hamilton's equations are only valid at some initial time for the Moyal evolution. This is a conse- 
quence of two related concepts, namely that Moyal evolution is neither a canonical nor a coordinate (local) 
transformation. A necessary condition for both evolutions to coincide is that Moyal (and classical) evo- 
lution be both a canonical and a unitary transformation. Whether this is also sufficient deserves further 
investigation. Indeed if this condition were sufficient for coincidence of the two evolutions, then in a first 
instance it would suffice to analyze the classical solution only in order to infer beforehand whether the 
Moyal evolution is equal to its classical counterpart or not. There are a few comments that can be made 
in this respect. Let us assume for a while that there exist Moyal trajectories which are simultaneously 
canonical transformations and distinct from the classical counterparts. Let us consider equations (40): 

Q = Fgiq,p,t), P = Fpiq,p,t), (111) 

where Fi{z,t) = U{t) • [^;j,-H^] * U{—t). The solutions are such that Q = Q{q,p,t) and P = P(q,p,t). We 
may invert these relations and get q = q{Q,P,t), p = p{Q,P,t). Substituting these equations in (111) we 
get: 

Q = GgiQ,P,t), P = GpiQ,P,t), (112) 

where Gi{Q,P,t) = Fi{q{Q, P,t),p{Q, P,t),t). In classical mechanics one would expect: Gg{Q,P,t) = 
Gq{Q,P,0) = Fq{Q,P,0) = -^H{Q,P), etc. If the Moyal evolution is a canonical transformation then 
there should exist some function He{Q,P) (an effective Hamiltonian) such that: 






G<iiQ^P) = ^, G,iQ,P) = -^. (113) 



Then the Moyal equations would be equivalent to: 

^ 9P ' dQ- ^ ' 

Since at t = to we have HeIQ, P,t = to) = H{q,p)^ we conclude that if He is not explicitly time dependent, 
then, in view of (114), Moyal and classical evolutions would have to coincide. At this point it is not clear 
to us whether there exist Moyal trajectories which are a canonical transformation and different from the 
classical trajectories generated by the same Hamiltonian. What can be said from the previous analysis 
is that if such trajectories do exist then they are equivalently generated by a classical flow with a time 
dependent Hamiltonian. 
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Another byproduct of our analysis is following conjecture: whenever a classical system can be solved exactly 
in the space of smooth phase space functions with the usual commutative product • and Poisson bracket 
{, }, then the corresponding Moyal flow can be "solved^ in the space of ^-functions with noncommutative 
^-product and Moyal bracket [, ] . One just has to deform the classical solution by introducing judiciously *- 
products (cf.(81,82) and (95,96)). However, it is still unclear what is the precise prescription for substituting 
standard products by •-products in this context. An even if one finds such prescription the final result in 
terms of ^-functions is obviously not explicit: it is expressed in terms of ^-products which may be difficult 
(or impossible) to evaluate. The obvious question is whether one gains anything by doing so. Indeed, 
one may just as well express the solution in terms of •-products as in eq.(27). However, by inspection of 
our example 1, we realize two important facts. First of all in our solution (95,96) one has to evaluate the 
noncommutative exponential (/>± (94) for a quadratic variable qp, which is well known [9], [15], whereas in 
the universal solution (27) one needs to compute the noncommutative exponential for a quartic variable 
H oc q'^p^ (i.e. the formalism of •-functions considerably simplified the resolution of the problem). Secondly, 
if one is not capable of evaluating all the •-products in the solution, then our prescription is more apt for 
performing the 7i-expansions order- by-order. Indeed, for the formal solution (27), when one computes the 
7i-expansion, one has to bear in mind the fact that the Weyl symbol of the evolution operator U{t) may 
be singular in fi. 

Another fact that our example 1 illustrates is that the classical and Moyal trajectories may diverge 
dramatically during a finite lapse of time and that the perturbative /l-expansions tend to obscure the 
global behavior of the trajectories. If we consider the ordinary linear differential equations (63,74) derived 
in section 3, the usual existence and uniqueness theorems apply. Typically the range of validity of the 
solutions found at each order will be less restrictive than that of the global solution. This is obviously 
connected with questions of convergence of the series, when one re-sums all the terms. We feel that the 
suggested method of solving the equations first in the noncommutative algebra U-^ (A) may be a better 
starting point for analyzing the global aspects of Moyal trajectories. 



Acknowledgments 

This work was partially supported by the grants POCTI/MAT/45306/2002 and POCTI/0208/2003 of the 
Portuguese Science Foundation. 



Note added 

After this work was completed we have learned about the paper [30]. Several of the results of that paper 
are coincident with ours, although the two approaches seem to be complementary. In addition, we present 
several examples. 
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